Absence of Overscreened Kondo Effect in Ferromagnetic Host 
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We study the low temperature behavior of a bound- 
ary magnetic impurity S' = 1/2 in an open ferromagnetic 
Takhatajian-Babujian spin-S 1 chain. For antiferromagnetic 
Kondo coupling, it is shown via Bethe ansatz solution that 
the impurity spin is always locked into the critical behavior 
of the bulk. At low temperature, a local composite of spin- 
(5—1 /2) forms near the impurity site and its contribution to 
specific heat is of simple power law Ts . The absence of over- 
screened Kondo effect is due to the large correlation length of 
host spins which is divergent near the quantum critical point. 
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Multichannel Kondo problem was originally proposed 
to study the magnetic impurity behavior in real metalsu. 
At low temperature for antiferromagnetic Kondo cou- 
pling, the impurity spin S' is completely compensated by 
the surrounding n-channel electrons when 2S' — n; while 
it is undercompensated for 25" > n or overcompensated 
for 2S' < n. The later situation (so-called overscreened 
Kondo effect) bears no resemblance to a non-interacting 
gas of electronsa and provides us a possible interpretation 
of noijL-Fermi liquid behavior observed in certain uraniun 
alloysB as well as in certain tunneling problems in two- 
level systems! As this effect is immediately destroyed 
by weak anisotropy, it has been a challenge to deduce its 
signature relevant to other anomalous behavior in a va- 
riety of exotic metals, especially when the systems show 
reduced effective dimensionality, or undergo the proxim- 
ity to quantum critical phase transitions!]. In this paper 
we will concentrate on multichannel Kondo problem in 
a typical quantum critical ferromagnet, with the number 
of channels n being represented by the host spin, n = 2S. 
Contrary to the problem in antiferromagnetic host chain, 
where overscreened Kondo effect appears when S > S' , 
we will find that due to the divergent correlation length 
near the critical point, the overscreened Kondo effect 
does not exist in ferromagnetic host. 

While the Kondo problems in paramagnetic or antifer- 
romagnetic hosts have been intensively investigated over 
three decades, they attract little attentions in ferromag- 
netic host. The main reason is the difficulty in approach- 
ing this problem, i.e., usual perturbation techniques as 
well as bosonization fail in the critical region because of 
the absence of the "large energy scale" , i.e., the Fermi en- 
ergy. Recent investigations have been shown that as the 
system closes to a quantum critical point, the impurity 



behavior depends strongly on the host properties and is 
non-universaE. Another reason is due to the fact that the 
existence of Kondo effect in ferromagnetic host remains 
controversial. At a first glance, the traditional Kondo 
screening seems possible only at temperature higher than 
the Curie point where host magnetization vanishes. Oth- 
erwise, a local impurity could not change the total mag- 
netization, opposite to what happened in dilute alloys. 
However, as far as the screening of a local moment is 
concerned, the issue should be how the impurity behav- 
ior changes as the temperature lowers down. The above 
picture does not necessary hold true even for the host 
goes away from the critical point where the long range 
order is destroyed by strong fluctuations. Yet there is no 
rigorous proof of the absence of overscreening in low tem- 
perature ferromagnet . In a earlier publication, Larkin 
and Mel'nikov first studied the Kondo effect in an almost 
ferromagnetic met all With traditional perturbation the- 
ory they shown that the impurity susceptibility is almost 
Curie type with logarithmic corrections at intermediately 
low temperature and the usual Kondo screening is not 
effective. More recently, the authors studied a similar 
problem in ID, i.e., a boundary Kondo impurity of spin 
S' in S = 1/2 ferromagnetfl. The model provides us a 
first exact soluble example for the Kondo problem near a 
quantum critical point. Via Bethe ansatz solution it has 
been shown that at sufficient low temperature, the impu- 
rity spin is completely screened by 2S" host spins. Thus a 
local composite forms near the impurity site and shows a 
simple power law dependence of specific heat on temper- 
ature. Therefore, these interest results strongly motivate 
us to know what will happen for a magnetic impurity in a 
multichannel ferromagnetic host at low temperature and 
if there exhibits the conventional overscreened Kondo ef- 
fect. 

In order to study the multichannel Kondo problem in 
ferromagnet via exact solution, we consider a magnetic 
impurity S' — 1/2 coupled to an open spin-S 1 chain which 
is integrable by Bethe ansatz. A well known integrable 
extension, of-, isotropic S — 1/2 spin chain to arbitrary 
spin S isETB 
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Hs — Jq >, Q2s{Sj ■ Sj+t), 



(1) 



where Q2s{x) is a polynomial of degree 2S of SU{2) in- 
variant quantities x = Sj •Sj+i- The general construction 
for Qs(x) is given in Rcf.[10,ll], one recovers the stan- 
dard Heisenberg model and Takhatajian-Babujian model 
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for 5=1/2 and 5 = 1 respectively. Here we follow the 
definition of H$ appeared in Ref. [11] so that the single- 
magnon dispersion coincides with that of the standard 
Heisenberg ferromagnet. Thus for Jo < 0, the model 
(Q) is qualitatively equivalcnt-tp the standard spin-5 fer- 
romagnetic Heisenberg chaintll. The boundary impurity 
spin is now coupled to the first host spin as 
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J s x ■ 5'. 



(2) 



The novelty of this construction is that the total HaroiJta- 
nian Hs + Hi mp is again integrable for any Jo and JE3E± 
The result is compared with the problem of a single im- 
purity in periodic spin chainsEJ, where the integrability 
requires a fixed Kondo coupling J and a fine turned im- 
purity interaction instead of simple form (|J) . As long as 
spin dynamics is concerned, the integrable spin-5 chain 
represents the multichannel host, because Kondo screen- 
ing of the impurity is equivalent to usual multichannel 
Kondo problem in antiferromagnetic case( Jo > 0)o the 
impurity spin is completely screened for 5' = 5; partially 
screened for S' > 5, with Schottky anomaly when an ex- 
ternal magnetic field H is applied; and overscreened for 
S' < 5, giving rise to quantum critical behavior. 

The model Hamiltonian is solved within the fr. ame .- 
work of Bethe ansatz approach to impurity problem! 13 !! 15 !. 
Eigenstates of the model are parameterized by a set of 
rapidities {X a }^ =1 (M being the number of down spins), 
with total magnetization S z — NS + 1/2 — M. Each ra- 
pidity represents a magnon of energy — J 5/(A 2 + 5 2 ), 
they are solutions of Bethe ansatz equations 
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where e r (x) 



x+ir/2 

x — ir/2 ' 



eim P (x) = ei + 2c(x)ei-2c(x) with 



c = ^/(5 + 1/2) 2 — J / J. In the following, we only con- 
sider the case Jo < 0, J > 0, so c > 5 + 1/2. In the 
thermodynamics limit the solutions are classified by n- 
strings or n-magnons, they are associated with functions 
Tj n (X) ( with real rapidity Ac3), which in thermal equi- 
librium for finite temperature T and field H satisfy an 
infinite set of non-linearly coupled integral equations: 

hiTM = -^rS „, 2S + Gm[(l + f?n _ 1 )(l + r)n+i)}, (4) 

with the asymptotic condition limn^oo inr )nW _ njj^ = 
2xo- Here, G is an integral operator with kernel g(X) = 
1/2 cosh(7rA). Notice that the free energy, expressed in 
terms of r] n (X), consists of three parts F — F^uik + 
jj{Fimp + F e dge), bulk part is always the dominate one, 
impurity and edge parts are both of order of I/AT, edge 
part is only due to the open boundary condition and is 
irrelevant to spin dynamics. We will discuss bulk and 
impurity parts in the following. 

Ferromagnetic host. The groundstate of ferromagnetic 
host has M = with finite magnetization S per site. The 



pure states are those with spin aligned along an arbitrary 
axis. With finite T and H the elementary excitations are 
n-magnons, bound states carrying n quanta of reduction 
of magnetization along the ordered axis. Each n-magnon 
has energy e„ = 2| Jo|& nj 2s(A) with crystal momentum 
K = 4 J2p=i arctan " +2 f+|~ 2p . In open boundary prob- 
lem, K is positive. The energy spectrum e n (K) is a con- 
tinuous function, increases monotonically as K varies in 
[0, TTmin(n, 25)]. The long- wavelength limit, i.e., K ~ 0, 
e n {K) oj \Jq\SK 2 /n, corresponds to |A| ~ oo. The ab- 
sence of gap for all n-magnons at boundaries of Bril- 
louin zones nmin(n, 25^ indicates the exact cancellation 
of Umklapp processes^. Correspondingly in thermody- 
namics limit, the driving term in Eq.(^) dominates when 
T — > 0, all rjn(X) functions are relevant. The situation is 
quite different to antiferromagnetic host (Jo > 0), where 
only 7725(A) is relevant as T — > 0. So the solutions of 
Eq.(|j) are more complicated for J < . Here two 
limiting cases are important: (a) Weak coupling limit: 
T — ► 00. The driving term disappears, there is only one 
parameter in equations, xq — H/T, the solutions are 
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sinh 2 [(n + l)xo] 
sinh 2 xq 
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(5) 



for all A and n, this is the free spin limit, (b) Strong 
coupling limit: T — > 0. In this case the driving term 
diverges, ln(l + rj n ) ~ lnry„, we obtain 



r) n (\) = exp[(nff + 7ra n , 2 s(A)| Jo\)/T}. 



(6) 



At zero temperature, rj n = 00 for all n, i.e., no magnons, 
the ferromagnetic groundstate with S z /N = 5 and 
E/N — H is reproduced. At low but finite temperature, 
r\ n (A) as function of A and n show crossovers between the 
two asymptotic solutions (||) and (|^). In principle, they 
can be solved numerically for fixed T and H. To study 
the critical behavior of -S, = 1/2 ferromagnetic Heisen- 
berg chain, SchlottmannO proposed an analytic method 
to solve these equations, based on an elegant but sim- 
ple correlation- length approximation, aiadrtiie result co- 
incides with the numerical one very wellt-j'tS. His idea is 
now developed for arbitrary 5 as follows. 

Let us first notice that for sufficient large |A|, the driv- 
ing term becomes negligible and one reproduces the weak 
coupling solution Eq.(§) for all n, up to C(e^ |A| | J \/T). 
Similarly, for sufficient large string index n, rj n (X) is again 
close to Eq.(||), because the effect of the driving term is 
also negligible. On the other hand, for small |A| and n, 
the driving term dominates, rj n (X) approach the strong 
coupling solution (|^). For intermediate A and n we have 
a crossover region. By equating two solutions @,@ 
for small A and H, we define a crossover scale n c (T) as 
25|Jo|n c = T ln[n c (n c + l)]/2, or solving it iteratively for 
T«|J |, 
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The meaning of n c (T) is clear: for n > n c , rj n is close 
to weak coupling solution while for n < n c it is close to 
strong coupling one. Therefore n c (T) serves as the corre- 
lation length of host, because it is the average number of 
the correlated spins in thermal equilibrium. Similarly, for 
each n-magnon, we define a crossover scale for rapidity, 



2S\J \ 



ln(n + 1) T 



(8) 



Obviously AJj (T) is the correlation length in momentum 
space: for |A| > A^(T) (long wave-length or low en- 
ergy), ?/ n (A) is close to weak coupling solution, while 
for |A| < A^(T) (and n < n c ) it is close to strong 
coupling one. Thus we adopt the following strategy to 
solve rjn(X), or more conveniently, to calculate Cn(A) 
I'l l + r? n (A)] — [7ra ni 2s(A)| Jo| + nH]/T: (i) As a zero or- 
der approximation, we assume rj n (X) = oo for n < n c 
and |A| < A£ , and elsewhere is given by Eq.(||). (ii) 
Substituting this approximation into the right hand side 
of Eq.(|J), we get the first order approximation for Cn(A). 
The result is 



Cn(A)« £{ln[l 



^^yl-^I^CA^-A) (9) 
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sinh(l + n c )xo 
sinh n c xo 



with B mn (X c m - A) = (f_£ + J™)A mn (X - X')d\'. 
The leading contribution to bulk free energy Ff,^ — 
— T J dXg(X)(2s{ty for small T and H/T is estimated as 



Ft 



bulk 



-1 K — \1/2t 3 / 2 



-0.42|J 



(2SH) 2 



1 lnln(^ 



2S|Jp|- 



T 2 ln 2 (^%M) 
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The zero field dependence of the free energy, as given 
above, is due to the n-magnon contributions from those 
within the cut-off, n c ; while the contributions beyond the 
cut-off is of higher order, ~ T 2 . So this part eventually 
approach the exact one as n c — * oo or T — > 0. The 
finite-field contribution of the free energy, which comes 
from all kinds of n— magnons, is a result of competition 
of both parts, and it is approximately proportional to 
the correlation length n Cl giving rise logarithmic correc- 
tions to the bulk susceptibility. Notice that by rescal- 
ing 25 J — > J , the bulk free energy is proportional 
to that of the standard spin— 1/2 ferromagnetic Heisen- 
berg chain, so one obtains C huik = 25C 1/2 ~ T 1 ' 2 and 

Xbulk = 25xi/ 2 ~T- 2 ln- 1 T. 

Boundary impurity. When the impurity is coupled 
to the ferromagnetic host via antiferromagnetic J > 0, 
Eq.(|) has a pure imaginary mode A = i(c — 1/2). This 
mode contributes a negative energy Ei mp = | Jo\S/[S 2 — 
(c — 1/2) 2 ] < to the system, indicating the formation 



of a bound state, i.e., a local composite made of the 
impurity and its neighboring spin S. One notices that 
A = i(c — 1/2) is the only possible, boundary mode for 
S' = 1/2 and no boundary stringu is allowed because 
A = i(c+ 1/2) is forbidden as shown in Eq.(3). 

Therefore, by taking into account the pure imaginary 
mode, the BAE for the bulk rapidities read 



2N ( \ A e2c-3(Aq) 

e2c-HA a j 
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=n n zi(k±\ ). (ii) 



When 2c is an integer, the impurity free energy is the 
difference of contribution from two effective "ghost" spins 
S+ = c - 1/2 and S_ = c - 3/2B 

Fimp = ~\\ dA 5 (A)[C 2c -i(A) - C 2c - 3 (A)]. (12) 

One finds that the leading contributions from the im- 
purity is the same as that of ^— spin(rescaling of Jo is 
used), but is always negative, i.e., CV mp = — C1/2 ~ T 1 ^ 2 , 
Ximp = — X1/2 ~ (T 2 lnT) _1 , indicating the freezing of 
some degrees of freedom of the bulk. Notice that c- 
dependence arises only in the next orders. When 2c is 
not an integer, the situation is somewhat more compli- 
cated. In this case, there are additional contributions 
coming from the non-zero residue {2c} = 2c — [2c] ([2c] 
is the integer part of 2c), which renormalize the effective 
Kondo temperature Tk ~ 1/ cos(7r{2c}/2) as well as the 
physical quantities. Explicit calculation shows that they 
only change the subleading behavior. 

The effect of the local composite can be analyzed in 
the same way. This local composite is made of the im- 
purity spin and the neighboring host spin. In contrast to 
the single channel problem where the local composite is 
always a spin singlet, the present one is not a spin singlet, 
but with a residual spin S — 1/2. Its feature is encoded 
in BAE as follows 
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(A Q ) = exp{^(A)}[] ]"[ ei(Aa±A/j). (13) 



represents the phase 



where 0(A) = -i m [e 2 s (A)f^g}] 
shift of a bulk spin wave scattering off the local com- 
posite. In the limit J — > 00, one finds a momentum 
dependent phase shift 0(A) =/= indicating no complete 
compensation of the impurity spin in multichannel fer- 
romagnetic host, (except for 5=1/2 where 0(A) = 0). 
Similarly, For 2c ^integer, we estimate the local compos- 
ite contribution to the free energy 



Fu 



T 



dA 5 (A)[2C 2S (A) + C 2c - 3 (A) - Cac-i(A)]. (14) 



One immediately obtains C loc = (S - l/2)C 1/2 ~ T 1/2 , 
Xioc - (S - l/2)xi/ 2 - (T 2 lnT)-\ So the leading or- 
der behavior of the local composite is similar to that 



3 



of (S — l/2)-spin, indicating that the impurity is nei- 
ther completely screened as in one channel problem nor 
overscreened as in antiferromagnetic host, but is always 
locked into the critical behavior of the bulk. 

Discussions and conclusions. In this paper, we studied 
a spin- 1/2 boundary impurity coupled to the ferromag- 
netic spin-S* chain. The results were obtained with the 
assumption that the crossover from strong coupling to 
weak one is abrupt. A smooth crossover region should not 
modify the low-T and small H/T dependences of the free 
energies obtained. Moreover, if the correlation lengthens 
n c and A° are scaled, the low-T dependence of our results 
remain unchanged, only the amplitudes are rescaled. The 
leading zero-field dependence of the bulk free energy is 
a rigorous result, even the amplitude.-agrees well with 
those obtained numerically for S =pJ/2lI3, or analytically 
for S > 1/2 via spin- wave theoryE3. However, since all 
n— strings contribute, the situation is more complicated 
for the susceptibility, logarithmic corrections in suscepti- 
bility arise even without Kondo impurity, indicating the 
existence of a marginal variable which does not exist in 
classical spin wave solution. 

The impurity effects are controlled by the bulk proper- 
ties and the impurity coupling. In the case of antiferro- 
magnetic host, see Ref. [12,13], a boundary impurity will 
shows the usual Kondo effect when the impurity spin S' 
is less than or equal to the host spin S, i.e., S' > S, 
otherwise, it shows a typical overscreened Kondo effect, 
i.e., S > S' . Depending on the Kondo coupling </, or c, 
there exhibits some novel quantum critical phenomenon 
in both cases. In the case of ferromagnetic host, the 
single channel problem (S = 1/2, S' > 1/2) shows a 
novel Kondo diamagnetic effect, the impurity spin is com- 
pletely screened by surrounding 2S' host spins for anti- 
ferromagnetic J > 0. Other magnons see a local compos- 
ite formed around the impurity site and interact with it 
weakly. 

The situation for multichannel ferromagnetic is dra- 
matically changed. Due to the presence of antifrromag- 
netic coupled impurity, there is always a pure imaginary 
mode in addition to the bulk n-magnons. The magneti- 
zation is suppressed due to the local composite which is 
a spin S — 1/2 object. But the impurity spin is not over- 
screened, in contrast to that in antiferromagnetic host. 
Rather, it is locked into the critical behavior of the bulk, 
with impurity specific heat showing a simple power law 
Ci m p ~ . The local composite now effectively coupled 
ferromagnetically to the other host spins via quantum 
fluctuation, showing the same power law specific heat 
Cioc ~ T?. The absence of multichannel Kondo effect is 
due to the large correlation length n c which is divergent 
when T — > 0. Though these results are limited to ID , 
similar phenomena may exist in higher dimensions. In 
multichannel ferromagnet, the antiferromagnetic Kondo 
coupling always leads to a local composite of a smaller 
spin, and leads to ferrimagnetic state when there is a 
finite density distribution of the impurities. 
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